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ABSTRACT 

A recently-proposed technique, called the dimensional expansion, uses the space-time 
dimension D as an expansion parameter to extract nonperturbative results in quantum 
field theory. Here we apply dimensional-expansion methods to examine the Ising limit of a 
self-interacting scalar field theory. We compute the first few coefficients in the dimensional 
expansion for 72™, the renormalized 2n-point Green's function at zero momentum, for 
n = 2, 3, 4, and 5. Because the exact results for 72 n are known at D = 1 we can compare 
the predictions of the dimensional expansion at this value of D. We find typical errors of 
less than 5%. The radius of convergence of the dimensional expansion for / j2n appears to 
be ^rj- As a function of the space-time dimension D, 72™ appears to rise monotonically 
with increasing D and we conjecture that it becomes infinite at D = We presume 

that for values of D greater than this critical value, 72 n vanishes identically because the 
corresponding (j) 2n scalar quantum field theory is free for D > —\ . 
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In a recent letter[l] we proposed a new technique called the dimensional expansion, 
which can be used to obtain nonperturbative results in quantum field theory. The dimen- 
sional series uses the space-time dimension D as an expansion parameter. The first term in 
such an expansion is easy to obtain because a quantum field theory can be solved in closed 
form in zero-dimensional space-time. An advantage of dimensional expansions is that some 
of the nontrivial aspects of the interaction already appear at D = 0. (Traditional perturba- 
tive methods yield only noninteractive results in leading order.) The obvious question is 
how one can obtain the coefficients of higher powers of D. A detailed explanation of how 
to do so is given in a subsequent paper [2]. 

Here we use the dimensional expansion to compute the first four 72 n , the renormalized 
2n-point Green's functions at zero external momentum, for a self-interacting scalar quan- 
tum field theory in the Ising limit. Specifically, we calculate 74 to fourth order in powers 
of .D, 76 to fifth order in powers of D, 73 to sixth order in powers of D, and 710 to seventh 
order in powers of D: 

74 = -L[l + (1.180 ±0.001)15 + (0.620 ±0.001)D 2 + (0.18 ±0.02)D 3 

+ (0.03 ± 0.02)L> 4 + . . .] , 
76 = JL [14. (2.20 ±0.02)1) + (2.30 ± 0.03)D 2 + (1.50 ± 0.03)D 3 + (0.55 ±0.04)D 4 

+ (0.12±0.04)D 5 + ...] , 

78 = JL [14. (3.0 ±0.1)1)+ (4.5 ±0.1)D 2 + (4.2 ± 0.1)D 3 + (2.6 ± 0.1)D 4 

+ (1.2±0.2) J D 5 + (0.6±0.2)L> 6 + ...] , 

710 = — [1 + (4.11 ±0.02)1) + (8.0± 0.1)D 2 + (10.0 + 0.3)D 3 + (8.0± 0.3)D 4 
90 

+ (4.5 ± 0.3)D 5 + (1.8 ± 0.3)D 6 + (0.7 ± 0.4)D 7 + . . .] . 

(1) 

To obtain these dimensional expansions we use the graphical methods described in 
Ref. 2. These graphical methods rely on lattice strong-coupling techniques that were 
developed and explained in an earlier series of papers[3,4,5,6]. For the Lagrangian 

£ = \mx)f + \m^(xf + \g^xf (2) 

the Ising limit [7,8,9] is defined as the limit in which the unrenormalized coupling constant 
g tends to infinity while the renormalized mass, M, is held fixed. The Ising limit is 
conveniently obtained by choosing m 2 (x—g. In the limit g^oo the theory asymptotically 
approaches a two-state system. The Green's functions of this system are universal in the 
sense that they are independent of the power of (j> in the self-interaction term in (2); g(p 2k 
gives the same results as gcj) 4 for all k>2. 

Lattice strong-coupling methods are especially well suited to obtain the dimensional 
expansion of Green's functions in quantum field theory because the lattice integral for 
each graph is a polynomial in powers of D. This property leads to an efficient organi- 
zation of the graphs that contribute to each order in the D-series. We achieve a high 
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order in the graphical expansion by eliminating all graphs except those that contribute to 
the coefficients in the dimensional expansion under consideration. We employ an inter- 
mediate renormalization scheme to calculate the renormalized mass M and dimensionless 
renormalized 2n-point scattering amplitudes 72™ at zero momentum. We perform mass 
renormalization of the scattering amplitudes by eliminating the bare mass m in 72 n in 
favor of the renormalized mass M. We then use Pade extrapolation methods to derive 
a sequence of approximants for each coefficient in the dimensional expansions in (1) for 
each of the scattering amplitudes 72 n in the continuum limit. We believe that each Pade 
sequence gives an accurate approximation to the true coefficient in the dimensional expan- 
sion for 72n because these dimensional series are in good numerical agreement with known 
exact results for 72™ [3]. The series in (1) are exact at D = 0. At D = l the exact results are 
74 = \-> 76 = \i 78 = jqi 7io = jq: and the results for (1) are 74 = 0.250 ± 0.003 (1% error), 
76 = 0.26 ± 0.01 (4% error), 7s = 0.31 ± 0.01 (4% error), and 710 = 0.42 ± 0.02 (5% error). 

We obtain the graphical rules for the lattice strong-coupling expansion by observing 
that in the limit of large g the kinetic term in the Lagrangian (2) can be viewed as a small 
perturbation. Therefore, the generating function 

Z[J] = Af J V<p(x) exp{ -fd D x[±[d<j)(x)} 2 + \m 2 (t){xf + \g<t){xf - J{x)(j>{x)] } (3) 

for the quantum field theory associated with the Lagrangian (2) can be rewritten as 

Z[J) = ex P {\jd D xd D yj^V-\x-y)j^) Z [J] , (4) 

where V~ 1 (x — y) = d 2 S D (x — y) and 

Z [J] = Af Jvcj) exp{- fd D x[±m 2 (j)(x) 2 + \g(j)(x) 4 - J(x)0(x)]} . (5) 

The factorization in (4) of the partition function leads to the strong-coupling lattice 
expansion. By introducing a D-dimensional hypercubic lattice with lattice spacing a we 
rewrite (5) as 

00 

Z [J] = jVjJ I dt exp{-^a D m 2 t 2 -^a D gt 4 + a D J,t} . (6) 

1 -00 

Next, we expand in powers of Ji and, to obtain the Ising limit, we set 

2 2 — D /f7\ 

m = —aga , (7) 
where a is a dimensionless parameter considered to be small in the strong-coupling limit: 

00 °? 

zo[j] = ArUT,7^S aDji)2n dttn ~ 1/2e A~r 9 ^ ■ (8) 
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In the limit g — > oo the integral in (8) is asymptotic to a n multiplied by a constant 
independent of n which we absorb into Af. Thus, we write Z [J] in (8) as 

oo 1 

Z [J] = Nexp{a D YX£— ^Jf n V2n\} , (9) 

where the vertices are V 2 = a 2 a, V 4 = -2a 4+D a 2 , V 6 = 16a 6+2D a 3 , V 8 = -272a 8+3D a 4 , 
Vio = 7936a 10+4i:) a 5 and so on. The propagator on the lattice can be written in vector 
notation as V~ x = a-- D_2 [(l) - 2D(0)]. This notation was introduced in Ref. 4 where 
this discrete form of the propagator was used to evaluate lattice integrals. The lattice 
strong-coupling expansion is organized by the number of free propagators V~ l (in contrast 
to weak-coupling expansions where the number of vertices and not the number of lines 
determines the order). 

To compute 72 n it is necessary to calculate the one-particle-irreducible 2n-point func- 
tions K-2n for n = 1, 2, 3, 4, and 5, in the strong-coupling expansion and to find their 
Fourier transforms A2 n in momentum space at zero external momentum. We must also 
compute Q(p 2 ) A-2(j> 2 ) \ p 2=o to obtain the wave- function renormalization constant defined 

by Z~ x = 1+ q^ 2 ) A2 1 | p 2=o- We define the scattering amplitudes 72 n as the dimensionless 
renormalized one-particle-irreducible vertices at zero external momentum 

7 2„ = f£(0,0,...,0)M J '(»- 1 )- 2 » , (10) 

where M is the renormalized mass defined as M 2 = r^(0, 0). There are simple rules giving 
r 2n in terms of A 2m , m < n, which are explained in Ref. 4: 

r 2 = A2" 1 , 
r 4 = -A4A- 4 , 

T 6 = - A 6 A 2 " 6 + ^^A 2 A 2 - 7 , 



\o = - A 10 A 2 10 + i^A 8 A 4 A 2 11 + ^^A 2 A 2 11 - J^_A 6 A^A 2 
: V 6 A 2 A- 12 + , JgU^ AUr 13 + i^A 4 A; 



!0 ! A A 2 A -12 , 10' A 4 A -13 , A 4 A -13 

2(3!) 2 4! AeA4A2 + 2 (2!) 2 (3!) 2 4 2 + W 4 2 



We use the wave-function renormalization constant Z to renormalize the one-particle- 
irreducible vertices in an intermediate renormalization scheme according to r|^(0, . . . , 0) = 
Z n T2n(0, • • • , 0). In order to mass renormalize the scattering amplitudes 72 n , we eliminate 
the bare mass m, which is related to a through (7), in favor of the renormalized mass M. 
To that end, we simply invert the relation obtained for the renormalized mass 

M 2 a 2 = a- 1 -2D + (2D-^)a + (4D 2 -^D+^-)a 3 

o 9 2584 68164 . - (10) 

+ 32 D 3 - 132 D 2 + — — D-— — )a 5 {ll > 

15 945 

+ (-2048 D 5 - 4096 D 4 - 1024 D s - 800 D 2 + 480 D - 64)a 6 + . . . 
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to expand a in terms of y = a 2 M 2 : 

a = y -2Dy 2 + (4D 2 + 2D - ?-)y 3 + (-8 D 3 - 12 D 2 + 4 D)y 4 + (16 D 4 + 48 D 3 

-4D 2 - 14D+ — )y 5 + (-32D 5 - 160 D 4 D 3 + 140 D 2 - 52 D)y 6 

5 3 

636 

+ (64 D 6 + 480 D 5 + 560 L> 4 - 720 D 3 + 20 D 2 + 272 D - —)y 7 + . . . . 

(13) 

We then substitute (13) for a in every to obtain 

y D/2 

74 = ^— [l + 4Dy + (4D 2 -10I% 2 + 16lV 5 + (-80I) 2 + 30I))y 4 + (256 J D 3 

J. 

+ 104 D 2 - 192 D)y 5 + (-704 D 4 - 1736 L> 3 + 2508 D 2 - 656 -D)y 6 
+ (1792 D 5 + 10432 D 4 - 11232 D 3 - 3872 D 2 + 4992 D)y 7 + . . .] , (14) 

76 = [l + 6Dy + (12D 2 -6D)y 2 + (8D 3 -12D 2 -20D)y 3 

+ (48 D 2 + 48 D)y 4 + (-96 D 3 - 816 D 2 + 528 D)y 5 
+ (192 D 4 + 4640 L> 3 - 2736 D 2 - 560 D)y 6 

+ (-384 D 5 - 18432 D 4 - 10800 D 3 + 46512 D 2 - 23040 D)y 7 + ...], (15) 

3D/2 

7 8 = ^— - [l + 8D ? / + (24D 2 -8D) ? / 2 + (32D 3 -32D 2 )|/ 3 
56 

+ (16 D 4 - 32 D 3 - 36 D 2 - 18 D)y 4 + (896 D 2 - 448 D)y 5 
+ (-4192 D 3 - 920 D 2 + 2816 D)y 6 

+ (13184 D 4 + 52064 D 3 -92800 D 2 + 38400 D)y 7 + . . .] , (16) 

V 2D 

7io = [1 + 10 Dy + (40D 2 -10D)y 2 + (80 D 3 -60 D 2 )y 3 

+ (80 D 4 - 120 D 3 + 30 D)y 4 + (32 D 5 - 80 D 4 - 300 D 2 + 108 L>)?/ 5 
+ (1280 D 3 + 5040 D 2 - 4240 D)y 6 

+ (—2560 D 4 — 64080 D 3 + 76880 D 2 — 23040 D)y 7 + . . .] . (17) 

The strong-coupling expansions in (14-17) were obtained by treating the dimensionless 
parameter a = —a D ~ 2 m 2 /g as small in the limit where the bare coupling g tends to infinity. 
The relation in (12) explicitly carries the assumption of smallness over to the parameter 
y. This justifies the reversion of (12) into (13) and the subsequent reexpansion of the 
scattering amplitudes 72™ in powers of y. Up to this point we have taken the lattice 
spacing a to be held fixed. We expect that in the continuum limit a — > our expressions 
for 72^ become the corresponding quantities of the continuum theory. The continuum 
limit is subtle because as a — > the parameter y that we have taken to be small actually 
becomes infinite. Hence, subsequent terms in this expansion for the scattering amplitudes 
72 n in (14-17) are increasingly singular as a series in powers of y in the limit where a — > 0. 
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We use Pade extrapolation techniques to extract information from perturbation se- 
ries like those in (14-17), where the perturbative parameter tends to infinity. The Pade 
extrapolation method employed here uses as input a perturbation series of the form 



f(y) = y r (c + c iy + c 2 y 2 + . . . ) (r ^ 0) 



(18) 



where we assume that f(oo) is finite. We first take the rth root of both sides of (18) and 
divide by y to obtain 



Then we take the iVth power of the right side of (19) for N = 1, 2, 3, ... , reexpand and form 
the (0, iV)-Pade approximant. By extracting the coefficient of y N in the denominator of the 
(0, iV)-Pade and raising it to the power —r/N, we create a sequence of Pade extrapolants 
for /(oo)[10]. We apply this method to the scattering amplitudes 72 n in (14-17) and 
expand each (0, iV)-Pade approximant for all 72 n in powers of D. For each n, we obtain a 
sequence in N of (0, iV)-Pade approximants for each coefficient in the D-series of 72 n - In 
Fig. 1 we plot the (0, iV)-Pade extrapolants for the first four coefficients in the dimensional 
expansion of 74 as functions of 1/N. We indicate the errors in the Pade determination 
of the coefficients in (1). We truncate the dimensional expansions for each 72 n after that 
coefficient for which the estimated error of the following coefficient becomes significant 
compared to its absolute size. 

Observe that the series in (1) all have positive coefficients and therefore each 7 2n is a 
monotonically rising function of D. Each of these functions is plotted in Fig. 2. For each 
n, 72n+2 is growing faster than 72 n for increasing D. We believe that the radius of conver- 
gence of the .D-series for 72 n is likely to be D = , the space-time dimension for which the 
coupling constant g of a g4> 2n theory becomes dimensionless and the theory becomes renor- 
malizable. Since we expect that for values of D> ^£3, 72n vanishes[ll,12,13,14,15,16], we 
assume that there is a singularity (possibly a natural boundary) in the complex-!) plane 



f(y) 1/r 



(new power series in y) . 



(19) 



y 



at L> = 



In 



n-1 - 
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FIGURE CAPTIONS 



Figure 1. Pade extrapolation for the first four coefficients, bi (i = f, . . . ,4), in the 
dimensional expansion of 74 = ^(1 + biD + b2D 2 + bsD 3 + b^D 4: + . . .). For 
each coefficient bi we plot the value of the (0, iV)-Pade (shown as cross) 
as a function otl/N for N = 1, . . . , 11. The continuum value of each bi is 
the extrapolation of the sequence to N = oo. In Eq. (1) we list the results 
of this procedure for 74, 75, 7s, and 710. 

Figure 2. Plot of 74, 7 6 , 73, and 710 in Eq. (1) as functions of D. For each n, 7 2n 
rises monotonically and 72n+2 is growing faster than 72 n for increasing D. 
We believe that the radius of convergence of the D-series for each 72 n in 
Eq. (1) isD = ^. 
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This figure "figl-l.png" is available in "png" format from: 
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